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From an accurate determination of the inter-quark potential, one can study the running coupling constant for 
a range of values and hence estimate the scale Aj^. Detailed results are presented for SU{'i) pure gauge theory 
from a study of a 36* lattice at f) = 6.5. 
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1. Introduction 



In the continuum the potential between static 
quarks is known perturbatively to two loops in 
terms of the scale Aj^jg. For 5f7(3) colour, the 
continuum force is given by [[l| 



dy_ 



4:a{R) 
3 i?2 



with the effective coupling a{R) given by 
1 

47r[6olog(i?AK)-2 + ih/bo) loglog(i?Afl)- 



(1) 



(2) 



where bo = ll/lGTr^ and bi = 102 5g/121 are the 
usual coefficients in the perturbative expression 
for the /3-function and, neglecting quark loops in 
the vacuum, A 

lattice regularisation scale A^ 



1.048Aj^. Note that the usual 
0.03471A^s. 
At large separation R, the potential behaves 
as KR where K is the string tension. Thus in 
principle knowledge of the potential V{R) serves 
to determine the dimensionless ratio \/K/A 
which relates the perturbative scale A to a non- 
perturbative observable such as the string tension 
K. This is the basis of the method |^ used for 
SU{2) which we extend here to SU{3) colour. 

For SU{3) as for SU{2), the bare lattice cou- 
pling proves to be a poor guide to physical be- 
haviour in that asymptotic scaling to two loops 
is not yet manifest. However, the weaker require- 
ment of scaling is well satisfied: the dimensionless 
ratios of physical quantities are found to be inde- 
pendent of p. For example in SU{2), the poten- 
tial V{R) scales [| over a range of lattice spac- 
ing of a factor of 4 (from f3 = 2.4 to 2.85). That 
scaling but not asymptotic scaling is valid implies 



that the bare coupling constant derived from /3 is 
inappropriate and that an effective coupling con- 
stant derived from some physical quantity is a 
better choice. This has been emphasized by Lep- 
age and Mackenzie Q. It is also the basis of the 
method proposed by Liischer et al. |^ to extract 
the running coupling constant. 

Here we choose to determine the running cou- 
pling constant from the interquark potential be- 
tween static quarks at small distance R. This 
quantity can be determined in a straightforward 
way from lattice simulation on large volume lat- 
tices. Although we require small R and hence 
large energy 1/R to make most precise con- 
tact with the perturbative expression, the lattice 
method implies the presence of lattice artefacts 
when R^i a. Thus we need to work on the largest 
spatial lattice available, consistent with avoiding 
finite size effects. We present results from a 36^ 
lattice at /3 = 6.5 to achieve this These results 
are compared with previous UKQCD data |^] 
from 24^ x 48 lattices at /3 = 6.2 to check scaling. 

2. Lattice potentials 

Previous studies of glueballs and potentials in 
quenched SU (3) have found evidence |^ for finite 
size effects when the spatial lattice size L is such 
that mL < 9 where m is the 0++ glueball mass. 
This corresponds to ^/KL < 2.6. Our previous 
study on a 24^ x 48 lattice at /3 = 6.2 cor- 
responded to \fKL = 3.8 which appears to be 
comfortably in the large volume region where fi- 
nite size effects are negligible. For our study at 
larger /3 we are constrained by the memory limita- 
tion to a 36'* lattice and we selected (i ~ 6.5 since 
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this again corresponds to the same spatial size in 
physical units, namely \J~KL = 3.8. Even though 
direct study of the potential seems to suggest ^ 
that it is relatively insensitive to finite size efltects, 
we prefer to adopt the cautious policy of conduct- 
ing our study at a fixed large physical volume. 

The details of our analysis of the potential 
between static sources on these 36'* lattices are 
given elsewhere Q . Here we summarise the main 
method and results. We conducted a relatively 
high statistics study of the potential at large R 
to be able to get an accurate value of the string 
tension K. From a fit to i?-values from 4a to 
24a we obtain Ka} — 0.0114(4) where a system- 
atic error from T-extrapolation is included. The 
values of the force are collected in table |l]. 

In order to explore the potential at small i?, 
we measure at off-axis as well as on-axis separa- 
tions. This will allow us to measure the extent to 
which rotational invariance is broken down to the 
cubic invariance on the lattice. Furthermore, we 
find that a simple parametrisation in terms of a 
lattice one gluon exchange is able to parametrise 
successfully our results. This is non-trivial since 
there is only one free parameter to accommodate 
the breaking of rotational invariance and yet a 
satisfactory fit is obtained to all the off-axis and 
on-axis data. The details of the fitting proce- 
dure are exactly those pioneered successfully for 

From this successful parametrisation of the de- 
parture from rotational invariance, we can con- 
struct a corrected potential by replacing the lat- 
tice gluon contribution in the fit by the full 
continuum gluon propagator. These results are 
shown in table |^. 

3. Running Coupling 

It is now straightforward to extract the running 
coupling constant by using 

) = 4^^^^ R, R, 

where the error in using a finite difference is here 
negligible. This is shown in table || and is plotted 
in the figure versus Ry/K where K is taken from 
the large R fit. The interpretation of a as defined 



Table 1 

The force AV/AR at 
rived from T-ratio 4/3. 



average separation R de- 



R/a AV/AR 



z. 


/I 1 /t o 

,4i4z 


U.066/(4j 


3. 


,4142 


0.0344(4) 


4. 


,2361 


0.0286(8) 


5. 


,0645 


0.0226(5) 


5. 


,8284 


0.0196(19) 


6. 


,1623 


0.0190(16) 


6, 


,7678 


0.0174(9) 


7, 


,6056 


0.0161(10) 


8, 


,2426 


0.0161(11) 


9, 


,0000 


0.0149(3) 


11, 


,0000 


0.0130(4) 


13, 


,0000 


0.0134(5) 


15, 


,0000 


0.0121(5) 


17, 


,0000 


0.0132(6) 


19, 


,0000 


0.0117(6) 


21, 


,0000 


0.0125(8) 


23, 


,0000 


0.0126(7) 



Table 2 

The force AV/ AR and lattice artefact corrected 
force AVc/ AR at average separation R. The run- 
ning coupling a{R) derived from the corrected 
force is shown as well. The second error shown 
on a is 10% of the lattice artefact correction. 



R/a 


AV/AR 


AVrJAR 


a[R) 


1.2071 


0.2067(7) 


0.1607 


0.170(1)(5) 


1.7071 


0.0750(7) 


0.0930 


0.197(1)(4) 


2.1180 


0.0959(19) 


0.0664 


0.223(6)(10) 


2.5322 


0.0541(5) 


0.0523 


0.248(2)(1) 


2.9142 


0.0263(40) 


0.0424 


0.270(26)(10) 


3.0811 


0.0471(39) 


0.0368 


0.262(28)(7) 


3.3839 


0.0391(12) 


0.0371 


0.317(10)(2) 


3.9241 


0.0292(5) 


0.0290 


0.333(6)(1) 
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above as an effective running coupling constant 
is only justified at small R where the perturba- 
tive expression dominates. Also shown are the 
two- loop perturbative results for a{R) for differ- 
ent values of A^. 



Figure 1. The effective running coupling con- 
stant a{R) obtained from the force between static 
quarks at separation R. The scale is set by the 
string tension K. Data at /3 — 6.5 are from ta- 
ble ^ (diamonds) and at /3 = 6.2 (triangles). The 
dotted error bars represent an estimate of the sys- 
tematic error due to lattice artefact correction as 
described in the text. The curves are the two-loop 
perturbative expression with a(6.5)A/j = 0.060 
(dotted) and 0.070 (continuous). 



The figure clearly shows a running coupling 
constant. Moreover the result is consistent with 
the expected perturbative dependence on R at 
small R. There are systematic errors, however. 
At larger R, the perturbative two-loop expression 
will not be an accurate estimate of the measured 
potentials, while at smaller R, the lattice artefact 
corrections are relatively big. Setting the scale us- 
ing y/K = 0.44 GeV implies l/a{P = 6.5) = 4.13 
GeV, so i? < 4a(6.5) corresponds to values of 
> 1 GeV. This _R- region is expected to be ad- 
equately described by perturbation theory. An- 
other indication that perturbation theory is ac- 
curate at such i?- values is that AVc/AR at small 



R is found to be very much greater than the non- 
perturbative value K at large R. 

Even though the lattice artefact correction of 
all 6 off-axis points by one parameter is very en- 
couraging, the only way to be certain that lattice 
artefacts are eliminated is by the comparison of 
different (3 values (with thus different R/a val- 
ues at the same physical R value). Now this test 
was satisfied in an SU{2) study [HJ|]. Even so 
we can check independently in SU{'3) and we use 
UKQCD data § from 24^ x 48 lattices at /3 = 6.2. 
We follow the same procedure as above (for more 
details see 0) and obtain for the string tension 
Ka^ = 0.0251(8). This corresponds to a ratio of 
lattice spacings a(6.2)/a(6.5) = 1.484(35) to be 
compared with the two-loop perturbative ratio of 
1.404. From the small R results we then obtain 
a corrected force and hence a running coupling. 
Setting the scale from the measured string ten- 
sions, we also show the (3 — 6.2 results for this 
effective running coupling in the figure. There is 
excellent agreement with the results from [3 = 6.5. 

The easiest way to describe the value of the 
running coupling constant a is in terms of a scale 
or A value with the understanding that we are 
only determining a for a range of energy scales 
1/R - namely 1 to 3 GeV. The final estimate of 
A is made from the figure, weighting smaller R 
more heavily since the perturbative expression is 
more accurate as a{R) becomes smaller. We ex- 
clude the lowest R point since the lattice artefact 
correction is for i? > a. Remembering that the 
systematic errors due to lattice artefact correc- 
tion are estimates only and since these system- 
atic errors are dominant, we do not attempt a fit 
but we can conclude that our results are consis- 
tent with values of A lying in the range shown by 
the two curves plotted. From the data at /? = 
6.5, these curves have a(6.5)Ai^=0.070 and 0.060. 
Using the value of the string tension from the fit, 
we get ^/K / Kl~ 49.6(3.8). Moreover, this value 
is consistent with the evaluation at both /3 = 6.5 
and 6.2. A rather similar analysis by Bali and 
Schilling also agrees with this result. 
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4. Conclusions 

Using the bare coupling g derived from [3 — 
&/g^ and the two-loop perturbative relationship 
a{g) in terms of the scale gave | |To|,[7| ,pt the 
following slowly decreasing values of Va/A^ = 
93.0(7) and 96.7(1.6)(2.6) at /3 = 6.0; 85.9(1.5) 
and 86.4(1.0)(1.9) at /3 = 6.2 and 82.3(8)(1.7) 
at /3 = 6.4. Our present analysis sX fi — 6.5 
yields y/K/kL = 80.0(1.4). Clearly, the /? ^ oo 
limit lies below these values. Moreover the statis- 
tically significant decrease of these values is ev- 
idence that two-loop perturbative scaling is not 
obtained in terms of the bare coupling. Our 
present method which does not rely on the bare 
coupling gives the scaling result which should be 
independent of 13. Our estimate is y/K / Kl = 
49.6(3.8). This is sufficiently far below the values 
extracted from the bare coupling to imply that 
asymptotic scaling to two-loop perturbation the- 
ory is not "just around the corner" but will only 
be satisfied accurately at larger /3-values than 
those currently accessible to lattice simulation. 

Using an empirical definition of an effective 
coupling constant in terms of the measured 
plaquette action, it is possible to extrapolate 
\fK/ K L to a = and the resulting estimates | |Tl| , 
1^ are in agreement with our value. This provides 
additional support for the viewpoint that the bare 
lattice coupling has bad perturbative behaviour 
but that perturbative descriptions of small dis- 
tance phenomena at existing (3 values are reliable 
if a physical coupling definition is used. 

Our result for the running coupling Q;i{(i?) given 
in the figure and table |^ can be read directly as 
a|^((7) with q = 1/R since these schemes are so 
close to each other. Since we obtain results con- 
sistent with the perturbative evolution, we can es- 
timate the continuum ratio Vk/Ajj^ = 1.72(13) 
for pure SU(3) theory. Setting the scale using 
VK = 0.44 GeV, then gives Aj^=256(20) MeV. 
These results are obtained for rather modest en- 
ergies { 1/R ~ 1-3 GeV ) but there is evidence 
from studies in SU{2) where higher energies have 
been reached that the method is stable as the 
energy scale is increased somewhat. From lattice 
results for ratios of other non-perturbative quan- 
tities (glueball masses, critical temperature, etc) 



to the string tension, one can then determine their 
value in terms of Aj^ as well. 

Even though the scales probed in this work are 
relatively small (i.e. only 3 GeV), the agreement 
with the perturbative evolution of the coupling 
constant implies that this is a reasonable way to 
determine the coupling constant in terms of non- 
perturbative physical quantities. We are able to 
determine Aj^ relatively accurately compared to 
experiment. Of course experiment has full QCD 
with dynamical light quarks included while pre- 
cise lattice simulation of full QCD is still a con- 
siderable challenge. 
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